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Abstract: We analytically study two-color QCD with an even number of flavors at high 
baryon density. This theory is free from the fermion sign problem. Chiral symmetry is 
broken spontaneously by the diquark condensate. Based on the symmetry breaking pattern 
we construct the low-energy effective Lagrangian for the Nambu-Goldstone bosons. We 
identify a new epsilon-regime at high baryon density in which the quark mass dependence 
of the partition function can be determined exactly. We also derive Leutwyler-Smilga-type 
spectral sum rules for the complex eigenvalues of the Dirac operator in terms of the fermion 
gap. Our results can in principle be tested in lattice QCD simulations. 
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1. Introduction 

Quantum chromo dynamics (QCD) exhibits various phases under extreme conditions: the 
quark-gluon plasma phase at high temperature T and the color superconducting phase 
||l|, 1^ at low T and intermediate or large quark chemical potential The quark-gluon 
plasma phase is being studied experimentally in ultrarelativistic heavy-ion collisions and 
theoretically in lattice QCD simulations. On the other hand, the physics of color supercon- 
ductivity at intermediate which is important for the interior of neutron stars, is not yet 
fully understood, although many possible phases have been proposed, e.g., the meson con- 
densed phase, the crystalline Fulde-Ferrell-Larkin-Ovchinikov phase, the gluon condensed 
phase, etc. [||. The case of asymptotically large /i is better understood because in this 
region rigorous weak-coupling calculations can be performed owing to asymptotic freedom, 
and the ground state has been shown to be the most symmetric, color-flavor-locked (CFL) 
phase [Q. 

Given the enormous theoretical and phenomenological interest in QCD at nonzero 
density, it is unfortunate that first-principle lattice simulations are extremely difficult be- 
cause of the fermion sign problem: The Dirac determinant becomes complex for nonzero 
/i, and conventional Monte Carlo methods fail. However, there is a category of theories, 
including two-color QCD with fundamental fermions and QCD with an arbitrary number 
of colors and adjoint fermions, that are free from the sign problem because of an additional 
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anti-unitary symmetry. Although these theories are quite different from QCD in certain 
respects (e.g., the pattern of chiral symmetry breaking is different from that of real QCD), 
they also exhibit dynamical phenomena such as diquark condensation that are expected to 
occur in real QCD, and thus we hope that they can provide us with some clues about the 
phase structure of real QCD. Recent years have seen substantial progress in the understand- 
ing of these QCD-like theories at nonzero temperature and density, both analytically and 
numerically |, |, |, g |, |, 0, |ll|, 0, 0, |l5|, 0, |l8|, 0, ^, |2|, |^ . 



In particular, two-color QCD has been studied intensively. Two-color QCD at small fx 
has been investigated by chiral perturbation theory in the mean-field approximation [^, |^, 
later generalized in various directions |9|, |l^, |l^ , and many intriguing phenomena have been 
revealed, such as the occurrence of a second-order phase transition at /i = where m-,^ 

is the pion mass. Also, vector mesons have been included in the chiral Lagrangian, which 
led to the prediction of vector meson condensation [12, 2^, 2£]. However, these approaches 



are based on an expansion in small /x and cannot be used to probe the system at large fi. 

On the other hand, at asymptotically large /x the formation of a BCS superfluid of 
diquark pairs is expected based on the existence of an attractive channel between quarks 
near the Fermi surface. The gap for quasiparticles has been calculated perturbatively in 
||3^, In this case, chiral symmetry is spontaneously broken not by the chiral condensate 
but by the diquark condensate, and the pattern of chiral symmetry breaking is quite 
different from that at /i ~ 0. 

In this paper we study two-color QCD at large /i analytically. We construct the 
chiral Lagrangian for the Nambu-Goldstone (NG) bosons associated with the spontaneous 
breaking of chiral symmetry for the case of an even number of flavors and determine the 
form of the mass term by symmetries.^ Furthermore we consider the system on a finite 
torus and identify the corresponding "e-regime" (or "microscopic domain") of the theory 
in which the finite-volume physics is completely dominated by the zero-momentum modes 
of the NG bosons. In this regime we can calculate the finite- volume partition function 
analytically as a function of the light quark masses. From this result we derive a set of 
novel spectral sum rules for inverse powers of the Dirac eigenvalues. These sum rules 



are a generalization of the well-known Leutwyler-Smilga sum rules at /i = |35|, but 
our analysis is more rigorous in the sense that the occurrence of spontaneous symmetry 
breaking is not assumed, but derived, from the microscopic dynamics of QCD based on 
asymptotic freedom. We also propose that a microscopic spectral density can be defined 
similarly to the /i = case, and that it should be a universal function depending only on 
the global symmetries of the system. We hope that our exact results could be an important 
stepping stone to a deeper understanding of two-color QCD (and maybe even real QCD) 
at nonzero /i. 

Recently, spectral sum rules for inverse Dirac eigenvalues have been derived in the 
CFL phase of three-color QCD by two of us [^], but due to the sign problem it is not 
straightforward to test them directly on the lattice. On the other hand, the new sum rules 
obtained in this paper for two-color QCD can in principle be tested on the lattice. The 



^See |34j for a similar construction in the CFL phase of three-color QCD. 
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lattice spacing has to be sufficiently small compared to l/fi to avoid lattice artifacts, but 
this seems to be mainly a technical problem and not a fundamental one. 

This paper is organized as follows. In section ^ we give an overview of the symmetry 
breaking patterns depending on the presence of the chiral condensate, the diquark con- 
densate, and/or the chemical potential (see also the discussion in 0]). In section ^ after 
reviewing the hierarchy of scales near the Fermi surface, we construct the low-energy ef- 
fective Lagrangian for the Nambu-Goldstone (NG) bosons at large // based on the pattern 
of chiral symmetry breaking caused by the diquark condensate. We also derive a mass 
formula connecting the mass of the NG bosons to the quark mass and the fermion gap. In 
section ^ we define the new e-regime and exactly determine the quark mass dependence of 
the partition function. We then derive spectral sum rules for the inverse eigenvalues of the 
Dirac operator at large /i. We conclude in section |5|. In the appendix we give a detailed 
microscopic derivation of the mass terms appearing in the low-energy effective Lagrangian. 

2. Symmetry breaking patterns 

Note that in this paper we discuss the continuum theory. Symmetries and their breaking 
pattern on the lattice depend on the lattice regularization adopted for the fermions and 
must be analyzed separately ||5|, |^ . 

We first define our notation, which is based on the notation used in 0. The fermionic 
part of the Lagrangian in Euclidean space reads ■ip{'D{fi) + M.)il) with the //-dependent 
Dirac operator 

= 7^D^ + 70/i (2.1) 

and the mass term 

= ^(l + 75)M + i(l-75)Mt. (2.2) 

Here, -0 is a short-hand notation for Nf flavors of two-color Dirac spinor fields transforming 
in the fundamental representation of SU(2)coior- The are hermitian 7-matrices, with 
75 = 7o7i7273- The covariant derivative Dy is an anti-hermitian operator so that the 
eigenvalues of ^yD^ are purely imaginary. M is the Nj x Nj quark mass matrix, and for 
real and degenerate quark masses the mass term simplifies to M. = ml ■ 

Because of the pseudo-reality of SU(2) we have V{pi)T2C^^ = T2C^^V{^)* , where C is 
the charge conjugation operator and T2 is a generator of SU(2)coior- Together with chiral 
symmetry, {75,P(/i)} = 0, this can be used to show that if A is one of the eigenvalues of 
P(/i), so are —A, A*, and —A*. While this observation is sufficient to see that the fermion 
determinant det(P(/i) +7W) is real, it does not necessarily imply that two-color QCD is free 
of the sign problem; the point is that A and A* are degenerate if A is real [^. Consequently, 
a sign problem can arise in two-color QCD with an odd number of fiavors. We also note 
that in the simultaneous presence of a quark chemical potential ^ and an isospin chemical 
potential a sign problem can arise in two-color QCD with an arbitrary number of flavors 
Q. Therefore, in this paper we only consider an even number of flavors and /U/ = so that 
the sign problem is absent. 
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For M = and /i = 0, the fermionic part of the Lagrangian is symmetric under 
\J{2Nf). For fj, ^ 0, this symmetry is broken exphcitly to \J{Nf)L x \]{Nf)ji. In general, 
the chiral symmetry of the Lagrangian is not reahzed in the ground state of the theory. 
Below we list the symmetries realized in various phases of two-color QCD with Nj flavors 
of Dirac fermions in the fundamental representation ^.^^ 

1. fi = 0, {i^^) = 0, (VV) = 0: SU(2iV/) [D {'L2Nj)a] 

2. fi = 0, ii^^) + 0, (VV) = 0: Sp(2iVj) [D (Z2)a] 

3. / 0, (V^V) = 0, (VV') = 0: SU(iV/)L X SU(iV/)K x U(1)b [d (Z2jv,)a] 

4. ^ 0, (V^V) / 0, (VV') = 0: SU(iVj)v' X U(1)b [d (Za)^] 

5. /i # 0, (VSV') = 0, (VV') / 0: Sp(iV/)L X Sp(iV/)^j [D x (^2)^?] 

6. ^ 0, (V^V) + 0, (V'V) / 0: Sp(iVy)v/ [D (Z2)b] 

In the first two lines, U(1)b is contained in SU(2A''j) and Sp(2A''j), respectively. {'^2Nf)A 
is the anomaly- free subgroup of the axial U(l)^ symmetry. {TL2Nf)A is already contained 
both in {Z]\f^.)L x U(1)b and in {Zn^.)^ x U(1)b. Our notation is such that Sp(2) ~ SU(2), 
therefore Sp(n) is defined only if n is even. In the last two lines, an even number of flavors is 
assumed. The fifth line is the case relevant for this paper, and the corresponding symmetry 
breaking pattern will be discussed in more detail in section ^ 

The color-singlet diquark field is a short-hand notation for 

= Babi^lfC^^PHi , (2.3) 

where C is the charge-conjugation matrix, a,b are color indices, and i,j are flavor indices, 
respectively. The A'^^ x Nf symplectic matrix / is defined as 



1= . n - (2-4) 




where 1 is the {Nf/2) x (iVj/2) unit matrix. (For Nf = 2, / is equal to —i(J2-) As 
long as s-wave condensation is assumed, there may in principle be not only a color- and 
flavor-antisymmetric condensate but also a color- and flavor-symmetric condensate. How- 
ever, the single-gluon exchange interaction indicates that the antisymmetric condensate is 
energetically favored [^] , and thus we will not consider the symmetric one in the following. 

Besides the scalar diquark condensate ip'^C^^tp in (p^), there is in principle also a 
pseudo-scalar diquark condensate ifFC'ip. Although the two are not distinguished by a 



single-gluon exchange interaction, it has been shown [38, |3^ that the instanton- induced 
interaction favors the former. Based on this observation we also neglect the latter conden- 
sate in the following.^ 

^There is no SU(2) global anomaly |3^] since we deal with Dirac fermions. 

^The symplectic group Sp(2n) — {g £ U(2n) \ g^Ig = /}, where / is defined in (2.4). The dimension of 



Sp(2n) is 2n^ + n, and its center is Z2 irrespective of n. 

*In Q it was rigorously shown by QCD inequalities, which hold even at nonzero chemical potential owing 
to the positivity of the Euclidean path-integral measure of two-color QCD, that for massless quarks the 
lightest meson is the 0"*" diquark, and thus a condensation (if it occurs) must occur in the channel C'^^ip. 
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3. Effective theory for dense two-color QCD 



We now proceed to construct the low-energy effective Lagrangian for two-color QCD at 
large /x. A similar construction for the CFL phase of three-color QCD has been performed 



in m, H, |3|, |4| 



3.1 Hierarchy of scales 

We first give a summary of the energy spectrum near the Fermi surface. For simplicity, 
the discussion in this subsection will mainly be restricted to the chiral limit (i.e., M = 0). 
According to perturbative calculations at large /x |30, ^] we have ~ ('0^) <C (V'^)-^ 



Therefore chiral symmetry is broken by the diquark condensate, not by the chiral conden- 
sate. The symmetry breaking pattern is (see section ^ 

S\J{Nf)L X S\J{Nf)R X U(1)b X U(1)a ^ Sp{Nf)L x Sp{Nf)n ■ (3.1) 



From the seminal renormalization group analysis near the Fermi surface |^^, the fermion 
gap A for SU(2)coior is found to be 
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A~/x5-^exp(^-— J , (3.2) 
where g = g{fj,) is the small running coupling obtained from the one-loop beta-function as 

''W'" (ll-/V,)ln(M/A,„,„) (3.3) 
We assume Nf < 11 to ensure asymptotic freedom. In particular, we have 

Asu{2) « A « /i (3.4) 

for sufficiently large /U. 

We first discuss the breaking of the U(1)b and U(l)^ symmetries. The U(1)b sym- 
metry is spontaneously broken, with an associated NG boson H that is gapless irre- 
spective of the current quark mass. The U(1)a symmetry is broken explicitly by the 
anomaly/instantons for small but this explicit symmetry breaking disappears as // — > oo 
because instantons are screened in this limit |]4l[| . In the latter case the U(1)a symmetry is 
broken spontaneously by the diquark condensate, and we have an associated NG boson rj' . 
At /i < oo the T]' acquires a mass from the anomaly, but for fi ^ Asu(2) the contribution of 
the anomaly to m^/ can be neglected, and the rj' can be treated as a pseudo-NG boson p2|| . 
In other words, at large it is a good approximation to treat U(l)yi as a non-anomalous 
symmetry and assume that the partition function is independent of the vacuum angle 9. 

We now turn to the breaking of SU(iV/)L x SU(iV/)ij to Sp(7V/)l x Sp(7V/)j?. The case 
of Nf = 2 is special since SU(2) ~ Sp(2). In this case there is no additional NG boson. 



' {ipip) — was observed in lattice simulations even at intermediate /i [|ll| 
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For Nf = 4 and larger, the spectrum contains additional NG bosons tt whose masses are 
lifted by a nonzero current quark mass. Their total number is 



2{N] - 1) - 2(iV|/2 + Nf/2) = N 



f 



Nf-2. 



(3.5) 



Are there other light particles? Quarks acquire a gap A and decouple from the dy- 
namics near the Fermi surface. What about SU(2) gluons? In two-color QCD, the diquark 
condensate (^) is a color singlet, so the SU(2) gluons do not acquire a mass through 
the Higgs mechanism. The absence below A of particles charged under SU(2) indicates 
that the medium is transparent for SU(2) gluons and that neither Debye screening nor 



Meissner effect occur |42|. For these reasons, the low-energy dynamics of gluons is simply 
described by the Lagrangian of SU(2) Yang-Mills theory. It has been found in [^] that 
the confinement scale of in-medium SU(2) gluodynamics is considerably diminished from 
the value ~ Agu(2) at /i = to 



SU(2) 



A exp 



2y27r n ^ 
11 



< A 



SU{2) 



(3.6) 



due to the polarization effect of the diquark condensate. Since n/ g/S. is large, see (^), it 
follows that gluons become almost gapless at asymptotically large /i. A more quantitative 



analysis |18] even suggests that gluons are lighter than the rj' for sufficiently large //. 
However, since their coupling to NG bosons does not show up at leading order of the QCD 
weak-coupling calculations |34, we assume in the following that gluons do not interact 
with NG bosons. 

The discussion of this section is summarized in table 0. The explanation of why = 



or oc m 



for large ^ will be given in section 3.2 





< ^SU(2) 


> ^SU(2) 


chiral symmetry 


spontaneously broken 


spontaneously broken 


2 

mi 


oc m 


or oc m? 


instantons 


abundant 


suppressed 


U(1)b symmetry 


intact 


spontaneously broken 


U(1)a symmetry 


anomalous 


spontaneously broken 


ry' 


heavy 


light 


condensate 


~ {ipTp) < Kv^v)! 


~ (v^V') < (V'V') 


mass gap (quarks) 


~ ^SU(2) 


~ A > Asu(2) 


mass gap (gluons) 


~ ^SU(2) 


~ ^SU(2) ^SU(2) 



Table 1: Qualitative differences in the physics of two-color QCD at large and small ^. m denotes 
the (degenerate) current quark mass. Nf is assumed to be even. See text for further explanations. 
In the second column, (V'V') l('0V')l only if A* ^ 



3.2 Chiral Lagrangian 

We now turn to the construction of the low-energy effective Lagrangian associated with 
the symmetry breaking pattern ( |3.1[ ). We introduce dimensionless color-singlet Nj x Nj 
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matrix fields Dl and Dji, where Nj is even and > 2 (the case Nf = 2 will be considered 
separately at the end of this subsection). In terms of the microscopic fields, they can be 
expressed as 

{DLr ~ ii^lyci^i , {Dnr ~ {^iyc^'r ■ (3.7) 

Under SU(A^/)l x SU(A''/)ii; x U(1)a x U(1)b the quarks transform as 

i^L ^ ^''^^^gii^L , i^R ^ e-'^'^-^^gnijR , (3.8) 

where gi E S\J{Nf)i [i = L,R) and the phases a and P are associated with the U(l)yi and 
U(1)b rotations, respectively. Thus 

Dl - 9lDl9W^''^'^ , Dr ^ gRDRgle-'^^'^-^^ . (3.9) 

It is convenient to split the fields into U(l) parts A, V and the rest, 

Dl = ^lA^V , Dr = J:rAV , (3.10) 

so that 

Ei^giJ^igf {i = L,R), A ^ Ae'^^" , V^Ve"'^. (3.11) 

The mass term tj^LMtpR + 'iIjrM'^iI)^ is invariant if we treat the mass matrix M as a spurion 
field and let 

M ^ SLM^lje^'" . (3.12) 

The effective Lagrangian composed of S^,, S/j, A, 1/, and M should be invariant under 
( |3lTD and (|3l2D . There is no invariant combination that contains an odd number of factors 



of M, in contrast to Q where the effective Lagrangian contained a term linear in M. This 
difference can be understood as follows. As Nf is even, we can take gi = —1 and gR = 1 
(or gi = \ and gR = —1), under which M transforms as M — > —M while S^,, S/j, A^ and 
V remain unchanged. Therefore M must appear in even powers. The point is that the 
chiral condensate {'tpL'4'R) is negligible at large /j.; otherwise a NG field S ~ 'tpi'^pR would 
appear that transforms under g^ = —1 as T, ^ —T, so that odd powers of M could appear 
in the Lagrangian in combination with S. (Strictly speaking, the assumption that Nj is 
even is not essential here. We can reach the same conclusion from the fact that the explicit 
breaking of U(1)a by instantons vanishes at large /i, as discussed in fS^ .) 
At O(M^) the real- valued invariant combination is uniquely found to be 

Tt{M^rM^J:[) + c.c. (3.13) 

No mass term appears for V reflecting the fact that \]{1)b is not violated by a nonzero 
quark mass. and T,r, which are decoupled from each other in the chiral limit, are now 
coupled by M. Replacing T^l^r by does not yield new invariants. This is because 

the diquark condensate is formed in a flavor antisymmetric channel, i.e., S^^^ = — S^^/j as 
already noted in section ||. In addition, we have SS^ = —EE* = Iat^ and det S = 1 for 
S G S\J{Nf)/Sp{Nf), which also greatly reduces the number of nontrivial invariants. 
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We conclude that the effective Lagrangian in Minkowski space-time, to lowest order in 
the derivatives and the quark masses, and assuming Nf > 2 and even, is given by^ 

r- = ^{\doV\' - vj,m\'} + ^{\doAf - v'^,\d.Af} (3.14) 

+ ^ Tr {|aoSL|2 - vl\di^L\^ + (L ^ i?)} - cA2|a2 Tt{MErM^^I) + c.c.} . 

Here, fn, fri'i and are the decay constants of H, t]', and vr, respectively, and the v's 
are the corresponding velocities originating from the absence of Lorentz invariance in the 
medium. The coefficient c is calculable using the method of [B4|] and found to be 



c = A . (3.15) 

The detailed derivation of (|3.15|) is given in appendix Note that £. has no dependence on 
the ^-parameter. The effective chemical potential of O(M^) (the so-called Bedaque-Schafer 
term in the CFL phase of three-color QCD [^^) is not displayed here since it is suppressed 
by ~ l//i at large fi. Also, its contribution to the finite- volume partition function starts at 
0{M'^), which is sufficiently small compared to the leading 0{M'^) term (see section ^). 
From ( p^ we can derive a mass formula for vr in the flavor-symmetric case M = 



mlj^j. For this purpose, we define the NG fields vr? (i = L, i?, a = 1, . . . , {N'j — Nf)/2 — 1) 
corresponding to the coset space '&V{N f)i/Si>{N f)i as^ 

= UilUf , Ui = exp ( '""l^ I , (3.16) 



where I is defined in (|2.4D and the X" are the generators of the coset SU(A^j)/Sp(A^j) 
satisfying the normalization Tt{X°' X^) = NfS""^ and the relation = I{X°')'^ . Since 
the mass term in ( ^.14| ) mixes ttl and ttr, it is necessary to diagonalize the mass matrix 
for 7r£ ^ to obtain the genuine mass eigenvalues by setting 

= ^{ttI + n%) , = - n%) . (3.17) 

The resulting mass formula for 11'* and 11" reads 

mn- = , /^m|„ = AcA^m^ . (3.18) 

The massless modes correspond to simultaneous rotations of vrj and tt^ in the same 
direction, while the massive modes 11° correspond to rotations in the opposite direction (a 
similar discussion can be found in |4^]).* For the massive modes H", ( 3.18| ) is of exactly the 



same form as the expression derived in [b2], except for the numerical factor of c. However, 



®The Wess-Zumino-Witten term is necessary for completeness of the theory, but it is irrelevant to the 
ensuing analysis and neglected in the following. 

'^For Nf = 2, UIU^ — (det U)I = / is constant, in accordance with the isomorphism SU(2) ~ Sp(2). 
^If the quark masses are nondegenerate, the If" modes acquire masses. 
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our "pions" for two colors are two-quark (qq) states, while those for three colors are four- 
quark (qqqq) states in |32]. 

Similarly, we can also derive a mass formula for the rj' boson associated with the 
spontaneous breaking of U(l)yi- We again assume the flavor-symmetric case M = mlNf 
The rj' field is defined by 

^ = exp(i^X — ). (3.19) 



Expanding to second order in terms of the rj' field in ( p. 14 ) yields 

/2 = AcA'^m'^ . (3.20) 



Since /^^^' ~ ^ [32| we obtain from ( |3.18| ) and ( 3.20 ) the relation 



mA 

mu,r,' '■ • (3.21) 

Finally, we consider the simpler case Nf = 2. In this case, the "pions" disappear. The 
first line of ( 3.14| ) requires no change while the second line is replaced by 



-c'A2|(detM)A2 + C.C.} with c' = ^ . (3.22) 



The above value of c' corrects the value of A/Sir'^ given in |jT^.^ 

4. Finite-volume analysis 
4.1 Microscopic domain 

In this subsection we identify the microscopic domain of dense two-color QCD and discuss 
some related intricacies. In the following, we neglect the H boson and the gluonic sector 
since they decouple from the dynamics of the NG bosons we are interested in and do not 
contribute to the quark mass dependence of the partition function. 

Let us study two-color QCD on a Euclidean four-dimensional torus P x L x L x L 
with linear spatial extent L and near-zero temperature /? = l/T ~ L. As is well known, 
the chiral limit and the thermodynamic limit do not commute in general. A particularly 
interesting regime to consider is the "e-regime" (or "microscopic domain" ) of the theory in 
which the zero-momentum modes of the NG bosons dominate the partition function |p6| . 
In the present setting, this regime is defined by the inequalities 

^ « L « , (4.1) 

where rn^j g ^, is the mass of the NG bosons (11, 11, rj') at large /i. The first inequality, 1/A <C 
L, means that the contribution of heavy non-NG bosons to the partition function can be ne- 
glected. The second inequality means that if the box size is much smaller than the Compton 



^We thank Thomas Schafer for a communication on this point. 
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wave length of the NG bosons, the functional integral over NG boson fields reduces to the 
zero-mode integral over the coset space S\J{Nf)L/Sp{Nf)L x S\J{Nf)ii/Sp{Nf)R x U(1)a- 
We briefly comment on the meaning of symmetry breaking assumed here. It is well 
known that a global symmetry cannot be broken spontaneously in a finite volume. How- 
ever, in two-color QCD the diquark condensate is a color singlet, which enables us to study 
diquark condensation phenomena by adding a gauge-invariant external symmetry-breaking 
field. This method is commonly used in lattice simulations, which are always restricted to 
a finite volume. This situation is in stark contrast to the CFL phase of three-color QCD, 
where the four-quark state composed of left- and right-handed diquark condensates is a 



color singlet, but the diquark condensate itself is not |l36 |. 
4.2 Massless spectral sum rules 

In this subsection we determine the dependence of the finite-volume partition function 
on the light quark masses and use it to derive a set of exact spectral sum rules for the 
eigenvalues of the Dirac operator 'D{^) defined in ( |2.lD . Let us first consider A'^^ > 4 (the 
case Nf = 2 will be discussed at the end of this subsection). Since the symmetry breaking 
pattern SU(n) Sp(n) is the same as at ^ = we can directly borrow calculational tools 



from P9| . The finite- volume partition function in the microscopic domain (O) is given by 



Z{M)= j dA j dTsL j dT.R exp T^cA^j^^ Tr(MS/jM^s[) c.c.} 

U(1)a S\](N;)L/Sp(Nf)L S\J{Ni)R/Sv{Nf)R 

= j dA j dUi j dUR ex.p -VicA^\^A^Tv{MURlUlM^UlIul)+c.c^ 



U(1)a SV{Nj)l SU(Af/)H 



(4.2) 



where Z is normalized to 1 in the chiral limit (i.e., all group volumes are defined to be unity) 
and V4 = is the Euclidean space-time volume. Note that Z has no ^-dependence. In 
going from the first to the second line of (4^), we have used S = UIU'^ as in ( |3.16| ), but for 



computational convenience we have chosen U £ SU(A^/) rather than U G S\J{Nf)/Sp{Nf). 
The integral over U is independent of the additional degrees of freedom we have introduced 
because UIU^ = / for [/ G Sp(iV/). 

From the invariance property of the Haar measure, ( |4.2D can be expanded as 

Z{M) = l + C2 Tr (M+M) + c'-l\Tr Alf + cf^ Tr (mI'M)^ + 0{M^) . (4.3) 

However, all terms of 0{M^^~'^) {k = 1,2,...) vanish because of the U(1)a integral. In 
particular C2 = 0, so the first nontrivial coefficients are c^^ and c;^\ both proportional to 
(V4cA^)^. To figure out their proportionality constants, we use the formulas given in |49, 
(4.11)-(4.17)], 

j dU exp { Tr(YUIU'^) + c.c. } = 1 + ^ Tr(Y^Y) + 0(Y^) (4.4) 

U{2N) 



^"Although there are some exceptions such as large- A'^ transitions and dynamical supersymmetry 
breaking 
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for any 2N x 2N antisymmetric matrix Y. (/ is given in (2^) but here has dimension 
2N.) Comparing ( |4^ ) and (p^ suggests the substitutions U AUr, N Nf/2, and 
Y -V^cA'^M^UIIUIm (Y is antisymmetric). We find 

Z(M) = J (IUl (V^cA^)^ ^ ^ ^ Ti{M^UlIUIm* ■ M'^UIIuIm)^ + 0{M^) 

= 1 + (y4cA2)2.^_A_|(TrMtM)2 - Tr(M^M)2| + O(M^) . (4.5) 

We note in passing that the above combination (Tr M^M)^ — Tr(M^M)^ also appeared in 
our previous work for the CFL phase |^6|. The above expansion is to be matched later 
against the one derived directly from the fundamental QCD Lagrangian. 



As a special case, let us consider the case of a flavor-symmetric mass term M = ml 



N 



f 



with real m. Now one of the two S\]{Nf) integrals (which corresponds to the integral over 
the massless modes 11" in ( |3.18| )) trivially drops out, leaving the far simpler expression 

Z{m)= j dU exp -V^cA'^m'^i^TriUIU^ I) + c.c.^ . (4.6) 

UiNf) 



A small consistency check can be done using (4^) in ( |4.6| ) to obtain (4^) with M = mljv^. 
This expression is of the same form as the one at ;U = p9| , (4.11)], which can be rewritten 
as a Pfaffian (5.13)]. By the same token (with the replacements x — > AcA^Vi^m? , — > 0, 
and 2Nf Nf, respectively, in ||4^, (5.13)]) we readily obtain the simple expression^^ 

^(-)=(^^^Pf(A), (4.8) 

where A is an Nj x Nf antisymmetric matrix with entries 

Apg ^{q- p)Ip+,(4cAV4m2) , p,q= • • • , (4.9) 

and Ip+q denotes a modified Bessel function. The analytical calculation of (|4.2D for non- 
degenerate quark masses (as has been done elegantly at ^ = for Dyson index 13d = 2 



||50| , |5l| , p^ , |53| as well as for = 1, 4 |Q) is left for future work. 

Having determined the light quark mass dependence of the partition function, we can 
now derive spectral sum rules for two-color QCD at large ^. The philosophy is well-known: 
By adding a weak external field to a given physical system and measuring its response 
thereto, one can gain information about properties of the system in the absence of the 
external field. In the present situation, the external field is a non-zero quark mass, and we 
can learn something about the properties of the Dirac spectrum in the chiral limit. 



^^The Pfaffian of an n x n antisymmetric matrix X (for even n) is defined as 

^^^^^ " 2"/2(n/2)! ^ sgn(cr)X<,(i)<,(2) . . . X<,(„_i)<,(„) , (4.7) 

where the sum runs over all permutations cr of 1, 2, . . . ,n. Note that this definition differs by a factor of 
l/2"''2 from H, (A.8)]. 
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We start with the properly normahzed partition function expressed in terms of the 
microscopic degrees of freedom, 



Z{M)= I [VA]det\v{fj,) + ^-^M+- — 



e '^s 



[VA]det^fV{fj,)e-^^ , (4.10) 



where Sg = f d'^x F^^F^^/A: is the gluonic action. Let us denote the complex eigenvalues 
of T){ii) by iXn, where the are real for = 0. In accordance with the discussion in 
section 3^, we assume that at large /i there are no exact zero modes, i.e., all are 
nonzero. Since all eigenvalues come in pairs (zA^, — iA„), we have 

Z{M) = I [VA] n'det(A^ + M) e'^^ / / [VA] (H'^n)"^' e'^^ 



(4.11) 



where (and later denotes the product (sum) over all eigenvalues A„ with Re A„ > 0, 
and 

(O) ^ j [VA] O (n'A^)"^' / j [VA] {j^X^r)""'^-'^ ■ (4.12) 



Because of the pseudo-reality of SU(2), the measure in (4.12) is real and positive definite 
(for even Nf). Using the relation 



det(l + e) = 1 + Tre + -[(Tr e)^ _ Tr e^] + 0{e^) 



(4.13) 



Z{M) in (4.11) can be expanded in terms of the quark mass matrix M as 

>(Tr MtM)2 



Z(M) = l + (j;'^)TrMtM+(^' ^ 



\2 \2 



+ ^ (| J^'^ ^ { (Tr MtM)2 -Tr(MtM)2} +0(m6). (4.14) 



Comparing the formal expansion (|4.14| ) with (|4.5|), we obtain novel spectral sum rules for 
two-color QCD at large /i,^^ 



ri)^(E' ' 



E 



/ 1 
Ai 



A^ A^ 

m<n ™ " 



(4CAV4) 



y'- 



0. 



4(iV/ - 1)2 ' 



(4.15) 
(4.16) 



which can in principle be continued to arbitrary orders of M. Note again that we have 
assumed even Nj with Nf > 4. The sums appearing in ( 4.15| ) and ( 4.16| ) are in fact real- 
valued for any fixed gauge configuration, i.e., before averaging over gauge fields. This can 
be seen from the pseudo-reality of SU(2): If A is one of the eigenvalues, so is A*. 



^^Sum rules for complex Dirac eigenvalues have been obtained earlier in a different context, i.e., at small 
/i in the conventional microscopic domain m^jfj. <C 1/1/ <C Aqcd, see, e.g., [psl. 



- 12 - 



Let us now define the spectral density 



(4.17) 



wliich is positive definite since the path-integral measure is (and thus p allows for a nat- 
ural probabilistic interpretation). Note that the gauge average in (4.17) involves massless 
quarks, see ( 4.12[ ). The case of nonzero masses will be addressed in section 4.3. 
Then (|4l6|) becomes 



with Cj 



(4cA2l/4)2 J^^ A4 

{w eC \ Re -u; > 0} and 



4:{Nf - 1)2 



(4.18) 



vr 



TTZ 



3A2y4^V^3A^ 



(4.19) 



where the numerical factor 4c = 3/7r2 has been included in the definition to simplify some 
analytical results that will be derived in section |4l^ . Equation ( 4.18 ) implies the existence 
of a microscopic spectral density defined by 



Ps{z) = lim 



(4.20) 



Since the sum rule ( 4.1^ ) is determined solely by the symmetry breaking pattern induced 
by the diquark condensate, we expect that ( 4.20D is a universal function fully determined 
by the global symmetries and independent of the details of the microscopic interactions. 
Our conjecture is corroborated by various facts known at fi = |56]: 



• The Leutwyler-Smilga sum rules are satisfied to good accuracy by the Dirac eigen- 
values computed in the instanton liquid model, a phenomenological model of QCD. 

• The microscopic spectral correlation functions derived from chiral random matrix 
theory (xRMT) generate the Leutwyler-Smilga sum rules exactly. Furthermore, they 
coincide with those derived from the finite-volume partition function of QCD to lead- 
ing order in the e-regime, because xRMT and QCD have the same global symmetries. 

• The microscopic spectral correlation functions predicted by ^RMT have been verified 
in a number of lattice QCD simulations. 

It is thus natural to expect that the way the thermodynamic limit of the spectral density 
near zero is approached is universal at large p, too. Further insights would be obtained 
if we could construct the appropriate random matrix model, a topic we hope to address 
elsewhere. 

It is interesting to ask how the vanishing sum rules ( 4.15| ) can be understood in terms 
of the symmetries of the distribution of the eigenvalues A^- One might conjecture that the 
eigenvalues occur not only in quartets (A, —A, A*, —A*) but in octets (±A, ±iX, ±X*, ±iX*). 
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This would be a sufficient condition for all averages of the type A ^'^"''^ {k 
vanish. However, this idea is flawed, as can be understood from 



1,2, 



to 



1 



which follows from (|4.15| ) and ( [4.16| ). Thus the symmetries are not realized per configura- 
tion, but only after averaging. We therefore turn to the spectral density p{\). It satisfies 
p{X) = p{X*) 
p(A) = p{iX), 



p{—\) = p{—X*). We now show that ( 4.15 ) follows if p{\) also satisfies 



d'X 



p{iX) 



(,;,)4.-2- ,d-X^ = ^. (4.22) 

The assumption p{X) = p{iX) does not contradict (4.16). We speculate that it is indeed a 
property of p. In view of so many symmetries, it is tempting to assume that p{X) = p(|A|). 
Then 



d'X 



A4 



|A|d|A| 



p_m 



dip e' 



0. 



(4.23) 



which contradicts ( [4.16| ). Thus p{X) is anisotropic. 

We finally consider Nf = 2. According to ( p. 22 ), the microscopic partition function 
is given by 



Z{M)= j dA exp |y4c'A2 

U(1)a 

= /o(2^4c'A2|detM|) 



(detM)^^ + C.C.]} 



= 1 + y(Kic'A2)2|(Tr MtM)2 - Tr(M+M)2| + 0(M®) , (4.24) 

where we used the Cayley identity for 2x2 matrices. Comparing (|4.24| ) with ( 4.14 ) we are 
again led to the vanishing sum rules ( [4.15| ) and to 



(4.25) 



which again hints at the existence of the universal microscopic spectral density ( 4.20| ). Note 
that ( f4.25| ) is actually the same as ( [4.16| ) with A'^^ = 2 because c' = 2c. The latter equality 
follows from Tt{MIM^I) = -2detM for Nj = 2, see ([aIsI) and (1a!6|) . In other words, 
the analysis of this subsection for Nf > 4 holds for Nf = 2 as well if we set Ti^ = S/j = / 
(corresponding to the trivial coset space SU(2)/Sp(2)) and Ul = Ur = 1 in (|4.2|) , and 



C/^ = 1 in ( |4.5| ). Of course, the physics is different since there are no "pions" for Nj = 2. 

At first sight, the above results are counter-intuitive in the following sense. Since the 
iXn are the eigenvalues of the Dirac operator ^yDy + 70/i for asymptotically large p, one 
naively expects A^ ~ p, and the inverse moments of the eigenvalues should decrease for 
larger p. But our exact results imply exactly the opposite!^^ — What we have overlooked 
is the existence of the Fermi surface. The typical momentum is comparable to the Fermi 
momentum, hence ji,Du is not necessarily small compared to 70/u. Delicate cancellations 
between these two terms lead to the above nontrivial sum rules. 



A diverges as fi ^ 00, see (B.2). 
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4.3 Massive spectral sum rules 



It is known that the massless spectral sum rules at ;U = due to Leutwyler and Smilga 
can be generalized to the massive case by an appropriate rescaling of masses , and that 
the double-microscopic spectral correlations (called so because masses and eigenvalues are 
rescaled simultaneously as the volume is taken to infinity) derived from random matrix 
theory reproduce these sum rules |58, 59, 30 1. Here we show that the new spectral sum 
rules we derived at large can be generalized to the massive case in a similar fashion. Such 
an extension will be useful in future comparisons with lattice gauge simulations, where light 
dynamical fermions are computationally expensive. 

We begin with the simplest case, i.e., A'^^ = 2 with masses mi and m2- According to 
the result in the previous subsection, we have 

Z = Io(amim2) with a = 2c'A^V4 = ^ — 

= Jf [m]]X[{Xl + ml){Xl + ml)]e-'^, (4.26) 

n 

where is a mass-independent normalization factor. Below we represent expectation 
values with respect to the massive measure by ((• • • )). 

Differentiation of Z by the masses yields massive spectral sum rules, e.g.,^^ 

^ 1 '7 //Y^' 1 W {am2f h{amim2) - h{arnim2) , . 



d{m\) W^A^+mf// 4 lQ{amim2) 

1 d'^Z „ , 1 , 1 „ .,2 

Z 9(mf )a(m2) A2 + ^ A2 + ml II 4 

1 d'^Z _ iis-^i 1 

Z " {\l^m\){\l^m\) 

{a.m2f 3Io{amim2) — 4/2 (ami m2) + Ii{amim2) 



(4.29) 



96 Io(amim2) 
For degenerate masses (mi = m2 = m) we have Z = iQ^am^) so that 

Zd{m'^Y W ^ (A2 +m2)2^ A2 +m2; // 2 Io(am2) • V -^^J 

In the massless limit equations ( [4. 27 )^( 4. 30| ) are consistent with the massless sum rules 
( 4.15 ) and ( 4.25| ). In terms of rescaled dimensionless variables. 



A m. = m.^^ (4 31) 

TT vr 



*We repeatedly used the identities In{x) — [In-i{x) + In+i{x)]/2 and In{x)/x — [In^i{x) — In+i{x)]/2n. 
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equations (|4.27D -(^^ read 

E' 1 \\ rfln In(x) — l2(x) , . 
I \\ ^ ^ 3/0(3;) -4/2 (x) +J4(x) 

- E'r^^ + = '"^^.^^y'^ • (4-35) 

It is surprising that the right-hand side of ( |4.33| ) is independent of the masses. 
We now define the double-microscopic spectral density by 



(4.36) 



mi=mi— fixed 



As mentioned before, the rescaling is the same for Nf = 2 and Nf > 4 because c' = 2c, 

and thus ( 4.36| ) applies to all (even) Nf. From the above sum rules we can derive several 

(2) 

equalities characterizing ps , e.g., 

2 pf\z]m,m) ^ lo(m^) - hjrn^) . . 

^ (z2+m2)2 4/o(m2) ' ^^■'^'> 



which is obtained from ( 4.35 ) and from ( [4.33| ) in the limit mi, m2 m. The determination 

(Nf) 

of Ps itself requires further theoretical analysis and is left for future work. 

For larger Nf the explicit expressions become increasingly involved. For Nf = 4 with 
equal masses, we have 

^-^ zi -|- // 4z 

n " 

where 

Z = 3/o(y)2 - 4/1 (y)2 + 3/2(^)2 . (4.39) 

Analogously, 



E fz2 +m2^2 



(z^ -|- m2)2 2;2_|_jfj2y // 

(4.40) 

One can check that the massless sum rules ( [4.15| ) and (4.16) are reproduced correctly in the 
chiral limit. Many more sum rules can be derived from the partition function for different 
fermion masses, but we do not work them out explicitly. We only mention that the sum 
rules corresponding to different masses satisfy nontrivial consistency conditions as some of 
the masses are sent to infinity, owing to decoupling |]5^. 
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5. Conclusions 



In this paper we have constructed the low-energy effective Lagrangian for two-color QCD 
with an even number of flavors at large quark chemical potential /i based on the symmetry 
breaking pattern induced by the formation of a diquark condensate. For Nj = 2, there are 
two Nambu-Goldstone bosons H and r]' associated with the spontaneous breaking of the 
U(1)b and U(1)a symmetries, respectively. The explicit breaking of U(1)a vanishes at large 
/i due to the screening of instantons. For Nf > 4, there are additional NG bosons ("pions" 
ttl and VTij) associated with the chiral symmetry breaking pattern S\J{Nf)i x S\J{Nf)ji — > 
Sp{Nf)L X Sp{Nf)ji. We derived a mass formula for the NG bosons in the case of degenerate 
quark masses. From the mixing of ttl and ttr by the mass term in the effective Lagrangian 
we found that half of them are exactly massless and that the other half are massive, with 
a mass proportional to the quark mass as a consequence of the Z(2)/, x Z{2)ji symmetry 
of the diquark pairing. 

On the basis of the mass hierarchy near the Fermi surface, we have identified a new 
finite-volume e-regime for the superfluid phase at large where the mass scale of the non- 
NG modes is characterized by the fermion gap A. In this regime, we can exactly determine 
the quark mass dependence of the partition function from the effective theory. Matching 
this result against the two-color QCD partition function, we have derived novel spectral 
sum rules for inverse powers of the complex eigenvalues of the Dirac operator. Our sum 
rules explicitly show that the Dirac spectrum at large fi is governed by the fermion gap A, 
unlike the spectrum at low /u, which is governed by the chiral condensate as shown in |49|. 

A remarkable property of two-color QCD is that the fermion sign problem is absent 
at nonzero fi (for even Nf with pairwise degenerate quark masses) because of the pseudo- 
reality of SU(2). Therefore, our sum rules can in principle be checked in lattice QCD 
simulations. This is in contrast to real (three-color) QCD where the severe sign problem 
prevents us from observing the presumed color superconductivity, although similar spectral 
sum rules could be derived in the corresponding e-regime |3^. Testing our sum rules for two 
colors on the lattice would enable us not only to verify the existence of the BCS superfluid 
phase originating from diquark pairing, but also to determine the value of the BCS gap 
A at large fi numerically. This is important since previous studies of two-color QCD at 
nonzero fi have only been able to determine the magnitude of the diquark condensate, not 
the gap itself. 

One should note that the diquark condensate does not necessarily imply a nonzero 
BCS gap A. At lower fi, just above m,r/2, the diquark condensate is a strongly-coupled 
Bose-Einstein condensate (BEG) rather than the weakly-coupled BCS pairing expected at 
large [^], although there is the possibility of a crossover between these two states |9| (see, 
however, [^5[). Also, the diquark condensate and the BCS gap are intrinsically different in 
the context of the Dirac spectrum in that the former (the latter) is governed by the radial 
(the anisotropic two-dimensional) distribution of Dirac eigenvalues |61]. We are planning 
to embark on a lattice simulation of two-color QCD at large fj, to test the ideas proposed in 
this paper and to make a first step towards providing a direct signature of the BCS phase 
of two-color QCD (i.e., the gap around the Fermi surface) at large fi. 
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The BEC-BCS crossover, or the continuity between the hadronic phase and the BCS 



superfluid phase |62|, may be closely related to the formal similarity of the partition func- 
tions for degenerate quark masses at small /i, see (5.13)], and at large n, see ( ^^ ) 
obtained in this paper. In the case of real (three-color) QCD with three flavors, it was 
recently suggested that the axial anomaly (or the effect of instantons), which acts as an 
external field for the chiral condensate, can lead to a crossover between the hadronic phase 



and the color superconducting phase |63, E^, EM as an explicit realization of the quark- 



hadron continuity in |62]. However, for QCD with two colors and an even numbers of 
flavors, the axial anomaly never acts as an external field for the chiral condensate, hence 
another mechanism to account for the crossover phenomenon would be necessary. 

Finally, regarding the Dirac spectrum, it would be interesting to obtain an analogue 
of the Banks-Casher relation at large ^ and the concrete form of the microscopic spectral 
density in the e-regime we identified. In particular, it is a challenging problem to construct 
the corresponding random matrix model, which has turned out to be very successful at ^ = 
and small fi, to reproduce our spectral sum rules and to elucidate universal properties of 
the Dirac spectrum at large fi for both two and three colors. (Earlier random matrix model 
approaches at large fj, include |65, 36, 67, |l^, 68, |2^.) The partial quenching technique 



(see, e.g., ||6^, ^ |71|1) may well be relevant in this regard. A detailed analysis of these 



issues is deferred to future work. 
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A. Microscopic derivation of mass terms in the effective theory 

In this appendix we give the detailed derivation of ( |3.15| ) and ( p. 22] ). For this purpose, 
we calculate the shift of the vacuum energy due to the quark mass from the microscopic 
theory (QCD), which should be matched against the result obtained from the low-energy 
effective Lagrangian ( p. 14 ). The starting point is the mass term of the microscopic theory 
given in [0^] , 

/:mass = |^(<ic^;;i)(v^^,^c<^) [ir^riT^rM,kM,i\ + (l^r,m^m^), (a.i) 

where k, I (a, b, c, d) are flavor (color) indices and {A = 1, . . . , — 1) are SU(A'c) 
color generators {N^ = 2 here) with normalization Tr(T"^T^) = 6^^ /2. The diquark 
condensates (II and that are antisymmetric in flavor and color can be written as 

mLC^^iL) = e'^'hjdL , {^IrC^U) = e'^'h^dR , (A.2) 
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where / is defined in (\2.4\) . The expectation values of and dn, in the ground state can 
be evaluated in weak coupling. Generalizing the result of [^] to arbitrary Nc, we obtain 



6Nc /i^A 
Nc + l -Kg 



Using (A. 2), (A. 3), and the Fierz identity 



1 



1 



) ) = 2'^ ^ - 2iv: 



^ac^bd 



the shift of the vacuum energy due to the quark mass is given by 



3A^ 
47r2 



Tr(M/M^/) + (M ^ M'^) . 



For Nf = 2 we have lij = Sij, and (A. 5) reduces to 

3A2 



AS 



Nf=2 



27r2 



detM + (M ^ M^) 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



By comparing ( [A.5| ) or ( [A .61 ) with the vacuum energy obtained from the low-energy effective 
Lagrangian ( |3.14| ) (with T,l = = I) or (3.22), the values of the coefficients are found 
to be 3 3 

respectively, where c' corrects the value of 4/37r^ given in 
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